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Diurnal  and  semi-diurnal  Hough  components 
of  surface  pressure 


(1.  V.  Groves 
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England 


Abstract 

Diurnal  and  semi-diurnal  Hough  components  of  surface 
pressure  are  evaluated  by  classical  tidal  theory  for 
previously  presented  profiles  of  water  vapour  and  ozone 
heating.  Values  are  compared  with  the  observational 
results  of  Haurwitz  &  Cowley  (1973)  and  show  significant 
discrepancies  which  are  considered  to  indicate  the  need  for 
additional  heating  to  be  identified.  The  present 


calculations  are  in  satisfactory  agreement  with  those  for 
semi-diurnal  modes  evaluated  by  Walterscheid  et  al,  (1980) 
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1 .  Introduction 

Larly  attention  was  given  to  water  vapour  absorption 
of  solar  radiation  as  a  generating  source  of  atmospheric 
oscillations  by  Siebert  (1961)  and  likewise  for  ozone  by 
Butler  &  Small  (1963).  Based  on  these  two  sources  of 
heating,  detailed  evaluations  of  diurnal  and  semi-diurnal 
atmospheric  tidal  fields  and  surface  pressure  oscillations 
have  been  presented  by  Lindzen  (1967,  1968),  Chapman  &. 

Lindzen  (1970)  and  Lindzen  &  Hong  (197*0 • 

Using  a  newly-constructed  model  of  ozone  densities, 
vertical  profiles  of  the  diurnal  and  semi-diurnal  Hough 
modes  of  ozone  heating  have  been  calculated  (Groves,  1980a) 
which  differ  significantly  from  those  of  Chapman  &  Lindzen 
(197g)  and  related  papers.  In  a  re-evaluation  of  the 
Hough  components  of  water  vapour  heating  (Groves,  1980b) 
cloud-related  scattering  properties  have  been  introduced 
and  have  the  effect  of  nearly  doubling  the  heating  at  8km 
altitude  and  of  halving  it  near  the  surface.  In  the 
light  of  such  revisions  to  the  ozone  and  water  vapour 
heating,  the  need  is  apparent  for  a  corresponding  evaluation 
of  the  atmospheric  tidal  response. 

The  present  paper  deals  with  surface  pressure  oscillations 
generated  by  the  diurnal  and  semi-diurnal  components  of 
ozone  and  water  vapour  heating.  Calculations  are  based 


on  classical  tidal  theory  (Groves,  1981)  and  for  the 
semi-diurnal  component  are  compared  with  the  evaluations 
of  Walterscheid  et  al.  (1980)  which  are  based  on  both 
classical  and  non-classical  assumptions  and  on  revised 
water  vapour  and  ozone  heating  rates. 

One  of  the  objectives  of  this  paper  is  to  compare 
calculated  values  with  previously  reported  observational 
results.  The  earlier  calculations  of  Lindzen  (1967)  for 
the  diurnal  tide  were  found  to  considerably  underpredict 
the  empirical  formula  of  Haurwitz  (1969)  at  latitudes  of 
less  than  30°  (Chapman  &  Lindzen,  1970,  Pig.  3*26).  The 
analysis  of  barometric  data  by  Haurwitz  &  Cowley  (1973)  has 
provided  a  more  detailed  set  of  results  which  offer  the 
possibility  of  comparing  individual  diurnal  and  semi-diurnal 
Hough  components  of  surface  pressure  with  theoretical  values 
for  different  seasons.  Attention  will  be  confined  to 
modes  of  oscillation  that  are  independent  of  longitude,  i.e. 
to  migrating  modes,  as  it  is  for  these  that  observational 
results  are  available.  Diurnal  modes  are  designated  by 
(1,1, n),  or  more  concisely  (l,n);  where  n  =  1,  2,  3».».  is 
a  sequence  of  propagating  modes,  those  with  n  odd  being 
equatorially  symmetric,  and  n  =  -1,  -2,  -3,...  is  a  sequence 
of  trapped  modes,  those  with  n  even  being  equatorially 
symmetric.  Gemi-diurnal  modes  are  designated  by  (2,2,n), 


or  more  concisely  (2,n);  where  n  =  2,  3»  4,...  , 
n  even  being  equatorially  symmetric. 


those  with 
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2.  hethod  of  calculation 

The  surface  pressure  oscillation  of  mode  (m,n)  is 
expressed  in  the  notation  of  Groves  (1981)  as 

=  (Tw  ccn,«-l  9"  Tn  0^(0)  0) 

where  t'  is  local  time,  «“  is  the  angular  frequency  of  the 
oscillation  and  0^  is  the  corresponding  Hough  function  of 
colatitude  Results  are  presented  in  terms  of 

r*"  =  ?"*R  +  L  "P”*1  (a) 

which  is  calculated  from 

T“  =  !p  1^(0)  W 

where  Pj(0)  is  obtained  from  Equ.  10.14  of  Groves  (1981) 
and  =  P00ui0/(r  »  P00  being  surface  pressure  and  coQ 

the  Earth's  sidereal  rate  of  rotation.  Pressure  then 
attains  a  maximum  or  minimum  value  of  |P^I  0™  at  local 
time  tan-1(P^/P^R)  according  to  whether  0®  is  positive 
or  negative.  The  basic  atmospheric  profiles  involved 
in  this  calculation  are  those  of  the  Newtonian  cooling 
constant  which  will  be  taken  as  zero,  i.e.  all  forms  of 
dissipation  are  then  neglected,  and  the  pressure  scale 
height  for  which  the  values  adopted  at  0(2)26  km  are 
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a.  75,  fa .46,  8.14,  7.82,  7.55,  6.98,  6.63,  6.22,  3*90,  3.99, 

b.  14,  6.24,  6.39,  6.36  km  and  at  28(2)150  km  are  taken  from 
the  mean  CIRA  (1972).  The  Hough  functions  introduced  in 
(1)  are  fully  normalized  being  taken  positive  at  the  equator 
if  symmetric  and  increasing  with  latitude  at  the  equator  if 
anti-symmetric.  Results  are  presented  as  plots  of  P^ 
against  .  Values  for  the  mid-season  months  of  January, 
April,  July  and  October  are  numbered  1,  2,  3  and  4  respectively 
and  those  from  Haurwitz  &  Cowley  (1973),  which  are  for  the 
intervals  hay  to  August  (J),  November  to  February  (D)  and 

the  remaining  months  (E),  are  qualified  by  the  letters  in 
brackets.  Comparisons  may  therefore  be  made  on  an 
approximate  seasonal  basis  between  1  and  D,  2  or  4  and  E, 
and  3  and  J.  Values  are  tabulated  in  terms  of  amplitude 
and  phase  (local  time  of  maximum  value)  in  Tables  1  and  2. 


3.  Diurnal  modes 

For  diurnal  modes,  water  vapour  heating  provides  a  much 
greater  contribution  to  surface  pressure  than  ozone  heating 
(Table  2).  In  the  case  of  the  trapped  sequence  of  modes, 
this  is  a  consequence  of  the  trapping  of  the  ozone -generated 
oscillation  at  stratospheric-mesospheric  heights  where  air 
densities  are  lower  and  variations  have  relatively  little 
effect  on  surface  pressure.  In  the  case  of  propagating 


Season 


1 


2 


3 


4 


P 

Ph 

P 

Ph 

P 

Ph 

P 

Ph 

(1,1) 

A 

290 

5.1 

283 

5.1 

255 

4.9 

283 

5.1 

B 

231 

4.9 

245 

4.8 

208 

4.9 

241 

4.8 

(1,-1) 

A 

40 

9.9 

84 

8.0 

104 

6.7 

84 

8.0 

B 

- 

- 

- 

- 

- 

- 

- 

- 

(1,-2) 

A 

418 

5.1 

458 

5.3 

449 

5.5 

458 

5.3 

B 

160 

6.0 

164 

6.0 

153 

6.0 

167 

6.0 

(1,-4) 

A 

248 

5.0 

246 

5.1 

196 

4.9 

246 

5.1 

B 

79 

6.0 

82 

6.0 

75 

6.0 

81 

6.0 

(2,2) 

A 

1100 

9.6 

11 33 

9.7 

992 

9.8 

11 33 

9.7 

B 

728 

9.0 

745 

9.0 

682 

9.0 

765 

9.0 

C 

759 

9.0 

- 

- 

759 

9.0 

- 

- 

D 

893 

9.0 

- 

- 

894 

9.0 

- 

- 

(2,3) 

A 

59 

9.7 

65 

0.3 

87 

1.3 

65 

0.3 

B 

30 

2.2 

4 

8.0 

26 

8.2 

6 

2.0 

C 

30 

2.2 

- 

- 

38 

8.4 

- 

- 

D 

55 

8.8 

- 

- 

49 

2.9 

- 

- 

(2,4) 

A 

166 

3.3 

182 

3.6 

151 

3.6 

182 

3.6 

B 

22 

3.4 

28 

3.8 

22 

3.4 

30 

3.7 

C 

16 

4.0 

- 

- 

18 

3.8 

- 

- 

D 

48 

3.7 

- 

- 

51 

3.7 

- 

- 

Table  1. 

Diurnal 

(m=1 ) 

and 

semidiurnal  (m= 

=2)  Houp;h  modes 

of  surface 

amplitude.  P=|Pm| ,  and  phase,  Ph,  for  modes  (m.n).  A  are  from  Haurwitz  and 
Cowley  (1973);  B  are  the  present  calculations  from  Table  2  (total);  C  are  the 
classical  and  D  the  non-classical  evaluations  of  Walterscheid  et  al  (198o). 
Columns  1,2, 3. 4  refer  to  the  four  seasons.  In  A  they  are  designated  D,E,J,E 
and  in  C  and  D  seasons  1  and  3  are  designated  DJF  and  JJA  respectively.  In 
the  present  calculations.  B,  they  are  for  January.  April,  July  and  October. 

P  is  in  microbars  and  Ph  is  the  local  time  of  maximum  (Qa  >  0)  in  hours. 
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modes,  the  depth  of  the  heating  region  in  relation  to  the 
vertical  wavelength  of  oscillation  is  important;  that  of 
ozone  heating  being  sufficiently  great  to  give  rise  to  surface 
pressure  contributions  that  tend  to  cancel. 

3 .  1  Trapped  modes 

For  trapped  modes,  water  vapour  and  ozone  contributions 
to  surface  pressure  are  in  phase  and  their  combined  values 
are  plotted  in  Fig  1.  For  (1,-2)  there  is  a  discrepancy 
between  calculated  and  observed  amplitudes  of  a  factor  of 
about  3  and  phases  differ  slightly.  Such  a  large 
discrepancy  is  difficult  to  reconcile  with  inaccuracies  in 
either  observed  or  calculated  values  and  would  seem  to  imply 
taut  water  vapour  and  ozone  heating  define  only  a  part  of 
tne  tidal  generating  source. 

The  comparisons  between  observed  and  calculated  values 
in  Fig.  1  for  (1,-4)  show  the  same  general  picture  as  (1,-2), 
out  some  similarity  is  to  be  expected  as  a  consequence  of 
tne  analytical  technique  and  the  limited  data  at  high 
latitudes,  where  this  mode  predominates.  Of  still  less 
independent  significance  would  be  comparisons  for  the  (1,-6), 
(1,-0,  (1,-10)  modes  which  are  the  other  modes  of  this 
sequence  included  in  the  analysis  of  Haurwitz  &  Cowley  (1973). 

n  significant  asymmetry  between  the  N  and  S  hemispheres 
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rig.  1  Diurnal  modes.  P^  and  P^1  are  plotted  in  jib 

on  the  horizontal  and  vertical  axes  respectively 
lor  modes  (m,n),  m=l,  n=l ,-l ,-2,-4.  The  local 
time  o'  ^  tan  is  shown  on  each  axis. 

1,2, 3, ^  denote  values  for  the  four  seasons. 
Broken  lines  are  position  vectors  to  calculated 
values  for  January,  April,  July  and  October. 
Continuous  lines  are  position  vectors  to  the 
observationally  derived  values  D,  iS  or  J  taken 
from  haurwitz  and  Cowley  (1973)  • 
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is  apparent  in  the  results  ol'  Haurwitz  &  Cowley  (1973)  which 
is  represented  mainly  by  the  anti-symmetric  mode  (1,-1.)  whose 
values  are  shown  in  Fig.  i.  No  theoretical  values  are 
shown  ior  this  mode  as  water  vapour  and  ozone  heating 
generate  almost  negligible  amounts  of  less  than  2  jib.  It 
may  be  noted  that  observed  values  do  not  decrease  to  zero 
with  the  solar  declination  but  have  a  significant  equinoctial 
value,  E.  The  only  other  anti-symmetric  mode  analysed  by 
Haurwitz  it  Cowley  (1973)  is  the  (1,-3),  but  few  stations 
exist  at  latitudes  south  of  40°S  with  which  to  justify  its 
signif icance  as  an  independent  determination. 

3.2  Propagating  modes 

For  these  modes  water  vapour  and  ozone  contributions  to 
the  surface  pressure  differ  in  phase,  the  ozone  contribution 
being  a  much  smaller  one  as  shown  by  the  dashed  lines  in 
Fig.  1  for  the  (1,1)  mode.  In  contrast  to  the  comparisons 
made  in  0  3.1,  the  (1,1)  results  show  closer  agreement  between 
observation  and  theory:  calculated  amplitudes  are  not  far 
short  of  observed  values  and  there  is  close  agreement  in 
phase;  also  July  amplitudes  are  slightly  smaller  in  both 
cases.  The  (1,3),  although  analysed  by  Haurwitz  &  Cowley 
(1973),  is  not  included  here  as  the  mode  predominates  at 
low  latitude  and  it  is  doubtful  whether  the  density  of 
equatorial  stations  is  adequate  to  resolve  a  migrating 
oscillation  of  such  low-latitude  structure. 
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4.  Semi-diurnal  modes 

The  (2,2)  mode  of  surface  pressure  oscillation  commands 
special  attention  on  account  of  its  large  amplitude,  which  is 
close  to  1000  pb,  and  its  accuracy  of  better  than  100  ytb. 
following  the  analysis  of  Butler  &  Small  (1965),  ozone 
heating  has  been  recognized  as  its  main  source  of  generation. 
Water  vapour  generates  about  half  the  contribution  of  ozone 
heating  with  the  same  phase  as  shown  in  Fig.  2.  Discrepancies 
between  observation  and  theory  are  apparent  in  Fig.  2  in  both 
the  amplitude  and  phase  of  (2,2),  but  earlier  classical 
evaluations,  such  as  that  indicated  in  Fig.  2  by  L  (Lindzen, 
1966),  have  shown  agreement  in  amplitude  to  be  very  good 
and  it  has  been  the  discrepancy  in  phase  that  has  previously 
attracted  attention  (Lindzen  &  Blake,  1971;  Lindzen,  1978). 

The  reduced  (2,2)  amplitude  of  the  present  calculations  is 
supported  by  the  classical  evaluations  of  Walterscheid  et  al. 
(I960)  for  the  solstices  which  are  marked  by  in  Fig.  2. 
walterscheid  et  al.  (I960)  also  evaluate  surface  pressure 
modes  under  non-classical  assumptions,  account  being  taken 
of  mean  zonal  winds  and  meridional  temperature  gradients, 
and  obtain  the  increased  (2,2)  amplitude  marked  by  in 
Fig.  2.  This  value  still  underpredicts  the  observed 
solstitial  values  by  1  0%  and  20%  respectively,  and  notably 
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does  nothing  to  reduce  the  discrepancy  in  phase. 

For  the  (2,4)  mode,  calculated  values  fall  within  the 
filled-in  square  in  Fig.  2  for  all  four  seasons  and  are  in 
fairly  good  agreement  with  the  classically  derived  solstitial 
values  of  walterscheid  et  al.  (I960)  which  are  marked  by  Wc. 
Such  values  agree  with  the  phase  of  the  observed  E,  J  and 
D  values  but  have  only  a  fraction  of  their  amplitude. 

Larger  amplitudes  have  been  obtained  by  the  non-classical 
calculations  of  Walterscheid  et  al.  (1980),  as  shown  at 
in  Fig.  2,  but  these  are  still  only  about  a  third  of  those 
ooserved.  Values  for  (2,9)  and  (2,6)  have  been  presented 
oy  Haurwitz  &  Cowley  (1975) ♦  but  no  comparisons  with 
calculation  are  made  nere  on  account  of  the  limited  data 
available  at  high  latitudes.  Results  are  however  presented 
in  Fig.  2  for  the  leading  anti-symmetric  mode  (2,5)  at  the 
solstices  for  which  there  is  very  good  agreement  between 
the  January  values  calculated  here  and  that  of  Walterscheid 
et  al.  (I960),  denoted  respectively  by  1  and  1(Wq)  in  Fig.  2, 
and  likewise  for  the  July  values  denoted  respectively  by  5 
and  5 ( Wq ) •  These  pairs  of  values  are  in  opposite  quadrants 
corresponding  to  the  reversal  of  sign  of  the  solar  declination. 
The  observed  values  do  not  however  decrease  to  zero  with  solar 
declination  but  have  the  equinoctial  mean,  E.  Non-classical 
(2,5)  results  denoted  by  l(W^),  5(W^)  differ  appreciably  from 
the  classical  ones  l(vV^),  5(WC);  and  are  in  substantially 
better  agreement  with  observation. 
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5.  Discussion 

For  diurnal  and  semi-diurnal  frequencies  comparisons 
have  been  made  between  Hough  modes  of  surface  pressure 
calculated  for  water  vapour  and  ozone  heating  and  the 
observational  results  of  Haurwitz  &  Cowley  (1973)*  The 
accuracies  of  the  latter  suffer  to  some  unknown  extent  from 
the  uneven  distribution  of  observing  stations  which  are  mainly 
concentrated  in  land  areas  at  mid-latitudes.  Accordingly, 
comparisons  between  calculated  and  observed  values  need  to  be 
interpreted  with  caution,  and  comparisons  have  not  been 
included  in  this  paper  for  some  of  the  modes  treated  in  the 
Haurwitz  &  Cowley  (1973)  analysis  which  predominate  at  high 
latitudes  where  data  are  sparse.  The  determination  of 
migrating  diurnal  modes  is  particularly  affected  by  the 
station  distribution  as  local  and  migrating  parts  of  diurnal 
oscillations  are  of  comparable  magnitude,  whereas  for  the 
semi-diurnal  oscillation  this  is  not  the  case,  the  migrating 
part  being  dominant.  For  this  reason,  the  underprediction 
at  low  latitudes  of  Haurwitz 's  (1963)  empirical  formula  for 
the  diurnal  oscillation  (Chapman  &  Lindzen,  1970,  Fig.  3.26) 
was  previously  attributed  to  the  likely  inaccuracy  of 
determination  of  the  migrating  diurnal  oscillation  at  low 
latitudes  where  station  density  is  low  (Chapman  &  Lindzen, 
1970,  p.  170).  From  the  comparisons  of  the  present  paper 
in  terms  of  modes  (Fig.  1),  theory  and  observation  are  in 


quite  close  agreement  for  the  low-latitude  confined  (1,1) 
mode,  whereas  serious  underprediction  by  a  factor  of  nearly 
3  arises  with  the  (1,-2)  mode.  It  is  not  clear  why  the 
underprediction  should  be  associated  mainly  with  the  (1,-2). 

One  possibility  is  that  a  source  of  tidal  generation  additional 
to  the  adopted  water  vapour  and  ozone  heating  is  present  with 
a  latitudinal  distribution  that  contributes  substantially  to 
(1,-2)  and  relatively  little  to  (1,1).  If,  however,  the 
earlier  appeal  to  low  station  density  is  followed,  the 
limitations  would  have  to  be  at  high  latitudes  where  (1,-2) 
predominates  and  (1,1)  is  inconsequential  and  not  at  low  latitudes 
as  previously  suggested.  In  this  case,  some  2/3  of  the  observed 
(1,-2)  mode  would  need  to  be  attributed  to  inadequately  resolved 
non-migrating  modes,  i.e.  (l,s,-2)  with  s  taking  values  close 
to  1;  and  the  problem  of  identifying  an  additional  source  of 


tidal  generation  would  still  remain. 

For  the  (1,-1)  mode,  water  vapour  and  ozone  heating 
provide  negligible  contributions  to  surface  pressure  and 
again  the  possibility  of  an  additional  tidal  source  attracts 
attention.  In  the  case  of  this  mode,  which  is  anti-symmetric, 
the  limited  distribution  of  S  hemisphere  stations  causes 
concern  about  the  accuracy  with  which  migrating  and 
non-migrating  oscillations  can  be  resolved  so  that  what 
was  analysed  and  identified  as  (1,-1)  may  actually  be  a 
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lumped  set  of  modes  (l,s,-l),  s  taking  values  near  to  1. 
liven  so,  the  source  of  excitation  of  such  modes  still  needs 
to  be  identified. 

The  (2,2)  is  the  largest  tidal  mode  of  surface  pressure 
and,  on  account  of  the  global  regularity  of  the  semi-diurnal 
oscillation,  is  the  best  observationally  determined  mode. 

The  discrepancy  between  observation  and  theory,  amounting 

2R 

to  about  400  |A.b  in  P2  (Fig.  2),  must  therefore  be  considered 
a  serious  shortcoming  in  the  application  of  tidal  theory. 

The  problem  has  been  considered  at  some  length  by  Lindzen 
&  Blake  (1971)  who  showed  that  neither  dissipation,  surface 
heating  nor  mean  winds  could  provide  an  answer.  A  later 
analysis  (Lindzen,  1978)  supported  the  suggestion  that  the 
discrepancy  may  be  due  to  the  release  of  the  latent  heat  of 
a  semi-diurnal  oscillation  in  tropical  rainfall,  although 
the  origin  of  the  oscillation  in  rainfall  remained 
unaccounted  for.  In  view  of  the  magnitude  of  the 
discrepancy  in  (2,2),  a  discrepancy  in  (2,4)  of  the 
magnitude  shown  in  Fig.  2  is  not  entirely  unexpected. 

The  observational  (2,3)  results  have  a  similar 
characteristic  to  those  of  the  (1,-1)  in  that  the  equinoctial 
mode  does  not  vanish  as  would  happen  if  the  heating  asymmetry 
depended  solely  on  solar  declination.  Also  like  the  (1,-1), 


observational  data  from  the  S  hemisphere  may  be  so  sparse 
that  the  results  actually  represent  a  lumped  set  of  modes, 
being  in  this  case  (2,s,3),  where  s  takes  values  near  to  2. 

The  evaluations  of  Walterscheid  et  al.  (1980)  of 
semi-diurnal  modes  by  classical  theory  are  denoted  by  Wc  in 
Fig.  2  and  show  satisfactory  agreement  with  the  present 
calculations.  The  effect  of  introducing  non-classical 
assumptions  was  examined  by  Walterscheid  et  al.  (1980)  and 
the  results  obtained,  which  are  denoted  by  in  Fig.  2,  are 
closer  to  the  observational  ones,  the  improvement  being  quite 
substantial  for  (2,3).  Significant  discrepancies  still 
however  remain  between  observation  and  theory  as  discussed 
above  and  indicate  the  need  for  additional  heating  to  be 
identified. 
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Terdiurnal  Hough  components  of  surface  pressure 
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England 


Abstract 

Profiles  of  terdiurnal  heating  due  respectively  to 
ozone  and  water  vapour  absorption  of  solar  radiation  are 
calculated  for  the  (3,3,3)  to  (3,3,6)  Hough  modes  and 
corresponding  surface  pressure  oscillations  are  evaluated 
by  classical  tidal  theory.  Comparisons  are  made  with 
earlier  evaluations  and  with  observationally  derived  results. 
The  calculated  solstitial  (3,3,4)  mode,  which  is  the  mode  of 
largest  surface  pressure  amplitude,  shows  good  agreement  in 
phase  with  observation  but  underpredicts  observed  amplitudes 
by  about  36  per  cent  in  contrast  to  earlier  evaluations 
which  were  based  on  a  now  unacceptable  basic  temperature 
profile. 


-  24  - 


CONTENTS 


Introduction 

Terdiurnal  Hough  functions 

Ozone  heating 

Water  vapour  heating 

Surface  pressure  oscillations 

Discussion 

Conclusions 

Acknowledgements 


page 

25 

26 
29 

34 

40 

45 

47 

48 


References 


49 


-  25  - 


1 .  Introduction 

Hann  (1918)  first  describe^.  the  terdiurnal  component 
of  surface  pressure,  its  main  characteristics  being  that  it 
undergoes  a  180°  phase  change  between  summer  and  winter  and 
its  amplitude  tends  to  zero  at  the  equinoxes.  Following 
Siebert's  (1961)  investigation  of  the  tidal  effects  of  the 
insolation  of  water  vapour,  Butler  &  Small  (1963)  carried 
out  a  similar  investigation  for  ozone  and  concluded  that 
nearly  all  of  the  terdiurnal  variation  of  surface  pressure 
could  be  accounted  for  by  these  two  sources  of  heating. 

Since  then  the  terdiurnal  atmospheric  tide  has  received 
little  attention  although  the  temperature  profiles  used  in 
the  earlier  analyses  have  been  modified  in  the  light  of  new 
data  and  more  detailed  water  vapour  and  ozone  models  are 
now  available.  A  re-evaluation  of  the  terdiurnal  surface 
pressure  therefore  seems  to  be  in  order. 

The  present  paper  calculates  Hough  functions  for  the 
(3»3»3)  to  (3,5*6)  modes  and  presents  their  form  graphically; 
the  relevant  heating  components  for  ozone  and  water  vapour 
are  computed;  the  surface  pressure  oscillations  arising 
from  this  heating  are  calculated  and-  comparisons  are  made 
with  previous  predictions  and  observed  results.  The  work 
is  an  extension  of  that  on  diurnal  and  semidiurnal  heating 
(Groves,  1980a, b)  and  corresponding  surface  pressures 
(Groves,  1981). 


2.  Terdiurnal  Hough  functions 

The  calculation  of  Hough  functions  has  been  described 
previously  in  detail  (Groves,  1979).  The  eigenvectors  {ar"$ 
corresponding  to  the  eigenvalues  X  of  a  matrix  whose  elements 
depend  on  the  Hough  mode  are  found  and  the  Hough  functions 
(£)  are  obtained  from 


©  =  £  arTr-,(/0 

where  p-  =  cos (colatitude ) ,  s  is  the  number  of  wavelengths 

that  fit  a  circle  of  latitude  and  o  is  the  normalized 

r,s 

associated  Legendre  function. 

The  matrix  can,  for  ease  of  calculation,  be  divided  into 
two  symmetrical  matrices,  one  producing  eigenvectors 
(°t&  ,°»  ^+2’°,  •••)  and  the  other  (o*.,  ,0,  V»,°,  •••)•  These 
eigenvectors  are  then  normalized  by  dividing  by  (a*-*'<£f^'  •  • 
and  (<L^+-  +  "-s  respectively.  Values  of  equivalent  depths 

obtained  from  h  =  X/0. 011349  km  for  modes  (3,?»n), 
n  =  3»  *  .  8»  are  12.890,  7*662,  5*083,  3*624,  2.714  and 
2.109  km  respectively.  These  values  agree  with  those 
which  were  given  by  Siebert  (1961)  to  one  less  significant 
figure  for  n  =  3»*  *  *  7*  The  following  series  expansions 


are  obtained  for  migrating  modes,  i.e.  for  s  =  3s 


©*  =  0.9965  P3  3  -  0.0857  Ps  ,  +  0.3571x10"^  P 


5,3 


7,3 


0.8807x10  P9  5  +  ... 


0*  =  0.9907  P4j5  -  0.1355  P6>5  +  0.8401xl0"2  Pg^5  - 


0.3025xl0~5  P10^3  + 


=  0.0830  P3  3  +  0.9798  P5>3  -  0.1813  Pr^5  +  0.0148  P^  ^  - 


0.7077xl0-5  P,,  ,  +  ... 

11,5 


0 l =  0.1339  P4>3  +  0.9652  P6>5  -  0.2233  P8>5  + 

0.0226  P 


10,3 


Comparison  of  the  coefficients  with  Siebert  (1961)  has  to 
be  made  with  the  aid  of  the  relation 

■P”  =  [VC^O4]  ^>0 

=  [2/C2^  +  »)]^R lWv  ^  —  0 

since  Siebert  used  the  serainormalized  Schmidt  form  P111  . 

n 

The  two  sets  of  coefficients  are  found  to  be  in  agreement. 
Table  1  shows  the  comparison  for  the  (3*3,4-)  mode.  The 
forms  of  the  (3,3,3)  to  (3,3,6)  Hough  functions  are  shown  in 
Fig.  1.  Signs  are  chosen  so  that  the  symmetric  functions 
(3,3,3)  and  (3,3,5)  are  positive  at  the  equator  and  the 
antisymmetric  functions  increase  with  latitude  at  the  equator 
going  north.  (3,3,3)  is  the  only  wholly-positive  terdiurnal 
Hough  function. 
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Fig.  1.  Terdiurnal  Hough  functions;  above  for  (3,3,3)  and  (3,3,4), 
below  for  (3,3,5)  and  (3,3,6). 
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3.  Ozone  heating 

The  calculations  are  based,  on  methods  and  models 
described  previously  (droves,  1980a).  The  ozone  model 
originally  extended  to  +60°  latitude  and  79  tan  height  (Lucas, 
1978)  and  was  derived  from  ground,  rocket  and  satellite 
measurements.  Extensions  to  +90°  and  80  km  were  then  added 
using  provisional  values.  Surface  albedo  values  of  0.07, 
0.14,  0.21  or  0.79  were  assigned  to  10°  latitude  x  10° 
longitude  areas  on  the  basis  of  all  sea,  all  land,  half  sea 
and  half  land  or  snow  cover.  do  daily  or  seasonal  variations 
of  surface  albedo  were  included.  Cloud  albedoes  of  0.25, 
0.629  or  0.679  were  given  to  10°  latitude  x  10°  longitude 
areas  on  a  seasonal  basis.  These  models  were  also  used  in 

the  calculation  of  water  vapour  heating,  results  of  which 
are  given  in  0  4. 

Table  1.  Series  coefficients  for  the  (3.3.4)  Hough  function 


Eigenvector 

Siebert ' s 
coefficients 

a  corrected  to 
Siebert 's  form 

a4 

0.99074 

0.66667 

1 

1 

a6 

-0.13991 

0.59470 

-0.164 

-0.1644 

a8 

0.00840 

0.48907 

0.012 

0.0117 

aio 

-0.00030 

0.43644 

-0.0005 

-0.0005 

Fig.  2.  Height  profiles  of  Hough  components  of  ozone  heating  for 

two  modes.  Key:  . (3,3,3);  - (3,3,5).  Butler  and  Small  (1963) 

spring  equinox  values  are  denoted  by  •  for  (3,3,3). 
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Fig.  3.  Height  profile  of  the  (3,3,4)  Hough  component  of  ozone  heating 
Butler  and  Small  (1963)  summer  solstice  values  are  denoted  by  •. 


Fig.  4.  Height  profile  of  the  (3,3,6)  Hough  component  of  ozone  heating 
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Table  2.  Hough  components  of  ozone  heating  associated  with  modes 
(3,3,3)  to  (3,3,6). 
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Hough  components  of  ozone  heating  were  calculated  from 
Equ.  3.21  of  Groves  (1980a)  and  are  presented  in  Figs.  2  to  4 
and  Table  2  for  the  (3,3,3)  to  (3,3,6)  modes.  In  all  cases 
maximum  heating  occurs  around  50  km.  For  the  (3,3,3)  mode 
the  July  maximum  is  some  6.5  per  cent  less  than  January's, 
in  keeping  with  the  increased  solar  distance.  For  the 
(3,3,4)  mode  the  peak  July  heating  is  1.7  times  that  of 
April  being  mainly  dependent  on  solar  declination  which  is 
9°  for  the  April  calculation. 

For  purposes  of  comparison,  the  values  obtained  by 
Butler  &  Small  (1963)  are  shown  in  Figs.  2  and  3.  The 
Butler  &.  Small  (1963)  values  are  generally  lower 
than  the  present  values  for  (3,3,3),  but  their  points  for 
(3,3*4)  in  July  are  in  good  agreement  with  the  present  curves. 

4.  Water  vapour  heating 

Hough  components  of  water  vapour  heating  are  calculated 
by  Equ.  6.9  of  Groves  (1980b)  and  include  the  effect  of 
scattering  by  clouds.  The  cloud  and  surface  albedo  models 
are  those  used  in  the  ozone  heating  analysis,  and  the  humidity 
model  was  based  on  the  maps  of  Wewell  et  al.  (1972)  for 
January,  April,  July  and  October  at  1000,  700,  500  and  400  mb. 

The  water  vapour  heating  profiles  are  presented  in 
Figs.  5  to  7  and  Table  3*  All  of  the  maxima  occur  at  around 
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Fig.  6.  Height  profiles  of  Hough  components  of  water  vapour 
heating  for  two  modes.  Key:  . (3»3»3) i  -  (3»3»5)» 


Fig.  7.  Height  profile  of  the  (3,3,6)  Hough  component  of  water 
vapour  heating. 
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t>  to  9  km  and  heating  decreases  rapidly  with  height  above 
the  maxima.  The  profiles  are  similar  in  form  to  the  diurnal 
and  semidiurnal  ones  (Groves,  1980b)  but  are  on  a  much  smaller 
scale.  maximum  heating  rates  in  mu/ kg  are  shown  in  Table  4 
and  contrast  with  the  (2,2,2)  October  peak  of  6.7  mW/kg.  As 
a  further  comparison,  the  ozone  (3,3,3)  October  peak  is  about 
60  times  greater  than  its  water  vapour  counterpart,  but  the 
two  resulting  surface  pressure  oscillations  (§  3)  are  in  the 
ratio  of  only  2  to  1. 


Table  Maximum  water  vapour  heating  rates,  mW/kg. 

Negative  signs  indicate  a  reversal  of  phase. 


(3,3,3) 

(3,3,4) 

(3,3,3) 

(3,3,6) 

January 

0.78 

-1.09 

r*H 

• 

o 

1 

0.38 

April 

0.70 

0.61 

-0.03 

-0.28 

July 

0.74- 

1.10 

-0.21 

-0.33 

October 

0.63 

-0.43 

-0.03 

0.22 
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8 .  Surface  pressure  oscillations 

Surface  pressure  oscillations  P  have  been  calculated 

o 

for  ozone  and  water  vapour  heating  profiles  by  classical 
tidal  theory,  the  method  of  analysis  being  that  previously 
employed  for  diurnal  and  semidiurnal  components  (Groves,  1981). 
For  the  (m,m,n)  mode,  PQ  is  expressed  as 

p0 -  (p:fwt'+  &,) 


where  t1  is  local  time  and  «r  is  the  angular  frequency  of 
oscillation.  Results  are  presented  in  Fig.  8  as  plots  of 
P™1  against  and  are  numbered  1,2,3  and  4  for  January, 
April,  July  and  October  respectively.  The  local  time 
S’-1  tan-1  (P^/P^)  is  shown  on  each  axis  in  Fig.  8. 
Amplitudes  and  phases  are  given  in  Tables  3  and  6. 

Comparisons  are  made  with  the  observed  values  quoted 
in  butler  &  Small  (1963)  which  are  indicated  on  Fig.  8  by 
A  for  the  annual  mean  and  J  for  a  solstitial  value. 
Comparisons  are  also  made  with  the  earlier  theoretical 
results  of  Siebert  (1961)  for  water  vapour  heating  and 
butler  it  Small  (1963)  for  ozone  heating,  the  combined  effects 
being  denoted  by  S  in  Fig.  8. 

For  the  (3»3»3)  mode,  the  amplitudes  of  the  results 
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.1,3 


(3,3,5) 


J - 1  nh 


•-2,4 


50 


Fig.  8  P®R  and  P™1  are  plotted  in  pb  on  the  horizontal  and  vertical 
axes  respectively  for  modes  (3,3,n),  n=3,4»5»  Local  times  of  maximum 

(0®  >  0)  are  shown  on  each  axis.  1,2, 3, 4  denote  calculated  values  for 

January,  April,  July,  October.  *  denotes  a  value  plotted  with  the 
sign  reversed.  S  denotes  the  combined  calculations  of  Siebert  (1961) 
and  Butler  and  Small  (1963).  A  is  the  observed  annual  mean  and  J 
the  observed  January  or  (-)  July  values  quoted  by  Butler  and  Small  (1963)* 


Season  1 

i  Ph  P 


(3,3,3) 


A 

- 

- 

- 

B 

38 

2.00 

35 

C 

- 

- 

48 

(3,3,4) 

A 

179 

1.89 

179 

B 

117 

1.98 

113 

C 

174 

2.11 

174 

(3,3,3) 

A 

- 

- 

- 

B 

2 

1.9 

3 
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• 

2 

4 

Ph 

P 

Ph 

P  Ph 

55 

1.78 

(annual  mean) 

2.00 

42 

2.00 

41  2.00 

2.00 

50 

2.00 

(annual  mean) 

-1.89 

29 

1.78 

(annual  mean) 

-2.02 

64 

-2.02 

55  1.98 

-2.11 

5 

-2.13 

(equinox) 

_ 

28 

1.78 

(annual  mean) 

2.0 

8 

-1.8 

8  -1.8 

Table  5.  Terdiurnal  Hough  modes  of  surface  pressure.  Values  are  of 
P  =  IP*1-  iP?1  I  in  ub  where  p3Rf  p3I 

are  defined  by  (5.1)  and  of  Ph, 

n  n  1  '  n  n 

the  local  time  of  maximum  ((*/  >  0)  in  hours,  for  modes  (3,3,n),  n=3,4»5» 

n 

A  are  the  observational  results  quoted  by  Butler  &  Small  (1963);  B  are 
the  present  calculations  from  Table  6  (total);  C  are  the  combined 
calculated  values  of  Siebert  (1961)  and  Butler  &  Small  (1963)*  Columns 
1,2, 3, 4  refer  to  the  four  seasons.  In  A,  1  and  3  are  for  January  and 

July  and  in  C  are  for  the  winter  and  summer  solstices.  In  the  present 
calculations, B,  1,2, 3, 4  refer  to  January,  April,  July  and  October. 
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PR 

Water 

P1 

vapour 

P  Ph 

PR 

Ozone 

P1  P 

Ph 

PR 

Total 

P1  P 

Ph 

(3,3,3) 

Jan 

0 

15 

15 

2.00 

0 

23 

23 

2.00 

0 

38 

38 

2.00 

Apr 

0 

15 

15 

2.00 

0 

27 

27 

2.00 

0 

42 

42 

2.00 

Jul 

0 

13 

14 

2.00 

0 

22 

22 

2.00 

0 

35 

35 

2.00 

Oct 

0 

14 

14 

2.00 

0 

27 

27 

2.00 

0 

41 

41 

2.00 

Ann 

0 

37 

37 

2.00 

0 

13 

13 

1.98 

0 

50 

50 

2.00 

(3,3,4) 

Jan 

0 

26 

26 

2.01 

2 

91 

91 

1.97 

2 

117 

117 

1.98 

Apr 

0 

-14 

14 

-1.99 

-1 

-50 

50 

-2.03 

-1 

—64 

64 

-2.02 

Jul 

0 

-2? 

27 

-1.99 

-2 

-86 

86 

-2.03 

-2 

-113 

113 

-2.02 

Oct 

0 

11 

11 

2.01 

1 

44 

44 

1.97 

1 

55 

55 

1.98 

Equ 

0 

0 

0 

- 

-1 

-5 

5 

-2.13 

-1 

-5 

5 

-2.13 

Sol 

0 

45 

45 

2.00 

15 

128 

129 

2.15 

15 

173 

174 

2.11 

(3,3,5) 

Jan 

0 

1 

1 

2.0 

0 

2 

2 

1.9 

0 

2 

2 

1.9 

Apr 

0 

0 

0 

- 

1 

-8 

8 

-1.8 

1 

-8 

8 

-1.8 

Jul 

0 

1 

1 

2.0 

0 

2 

2 

1.9 

0 

3 

3 

2.0 

Oct 

0 

i 

1 

2.0 

1 

-9 

9 

-1.8 

1 

-8 

8 

-1.8 

Table  6.  Terdiurnal  Hough  modes  of  surface  pressure  evaluated  for  water 
vapour  heating  and  ozone  heating.  PR,  P*  are  the  quantities  P^R,  P^ 
defined  by  (5.1)  and  are  in  ^b.  P  and  Ph  are  defined  in  Table  5*  The 

columns  under  Total  are  the  sum  of  the  water  vapour  and  ozone 
contributions.  The  entries  under  Jan,  Apr,  Jul,  Oct  are  the  present 
calculations.  Those  under  Ann  (annual  mean),  Equ  (equinox),  Sol  (solstice, 
January)  are  from  Siebert(l96l)  for  water  vapour  and  Butler  &  Small  (1963) 


for  ozone 


-  44 


(1,2,5  and  4  in  Fig*  8)  are  about  20  per  cent  less  than 
those  previously  calculated  (S),  but  are  still  close  to  the 
observed  value  (A),  being  within  about  15  |*b.  The  relative 
contributions  of  the  two  heat  sources  in  the  present 
calculations  are  however  substantially  different  from  the 
earlier  ones,  and  now  amount  to  65  per  cent  from  ozone  instead 
of  25  per  cent. 

The  antisymmetric  (3»3</+)  mode  attains  its  largest  values 
at  the  solstices.  Earlier  calculation  of  solstitial  values, 
denoted  by  S  in  Fig.  8,  has  shown  about  three-quarters  of 
the  amplitude  to  be  generated  by  ozone  heating  and  that 
together  the  two  heat  sources  account  adequately  for  the 
observed  value  (J).  In  the  present  calculation  ozone 
heating  again  accounts  for  about  three-quarters  of  the  total 
amplitude,  but  this  amount,  denoted  by  1  or  3*  in  Fig.  8,  is 
now  only  64  per  cent  of  that  observed. 

Fig.  8  includes  a  plot  of  (3»3,5)  values  although 
amplitudes  are  probably  too  small  for  any  useful  comparison 
to  be  made  between  theory  and  observation. 


L 
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6.  Discussion 

Previous  work  on  the  diurnal  and  semidiurnal  components 
of  ozone  and  water  vapour  heating  and  their  contributions 
to  surface  pressure  has  been  extended  to  the  terdiurnal 
component.  Although  much  smaller  than  the  diurnal  and 
semidiurnal  components,  the  terdiurnal  is  dominated  at  the 
solstices  by  its  (3»3i4)  mode  with  an  amplitude  of  179  pb 
and  comparisons  between  theory  and  observation  can  reasonably 
be  made.  Previous  calculations  by  Siebert  (1961)  for 
terdiurnal  water  vapour  heating  and  Butler  &  Small  (1963) 
for  terdiurnal  ozone  heating  have  given  a  total  surface 
pressure  oscillation  in  very  close  agreement  with  observation, 
not  only  for  the  (3»3»4)  mode  but  also  for  the  smaller  (3,3»3) 
mode  of  amplitude  55  p-b  (Fig.  8). 

The  present  results  for  the  (3»3»4)  mode  agree  closely 
with  the  previous  ones  in  phase  but  differ  in  amplitude 
accounting  for  only  64  per  cent  of  the  observational  value 
(Fig.  6).  Although  ozone  generates  much  (about  &)  of  the 
surface  pressure  oscillation,  the  decrease  is  not  readily 
attributed  to  changes  in  the  ozone  heating  profile  as  the 
present  (3, 3 >4)  solstitial  profile  and  the  earlier  values 
of  Butler  &  Small  (1963)  are  in  good  agreement  (Fig.  3)» 
Attention  has  therefore  been  given  to  the  effect  of  the 
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basic  temperature  profile  on  the  calculated  surface  pressure, 
nutier  it  Small  (1963)  adopted  a  temperature  profile  from 
nurgatroyd  (1937)  up  to  100  km  and  extended  it  to  130  tom  with 
a  thermospheric  gradient  of  4  K/km.  The  30  km  maximum  had 
a  value  of  about  295  K.  and  the  8o  km  minimum  203  K. 

Compared  with  the  profile  used  in  the  present  work,  the 
50  km  temperature  was  higher  by  23  K,  the  mesopause 
temperature  was  higher  by  8  K  and  lower  in  altitude  by  about 
10  km  and  the  thermospheric  gradient  rose  at  a  lower  rate. 

The  present  scale  heights  at  0(2)26  km  are  8.75»  8.46,  8.14, 
7.82,  7-53,  6.98,  6.63,  6.22,  5.90,  5-99,  6.14,  6.24,  6.39, 

6.56  km  and  at  28(2)150  km  are  taken  from  the  mean  CIBA  (1972). 
When  the  (3,3,9)  mode  is  re-calculated  with  the  present 
heating  and  the  basic  temperature  profile  used  by  Butler  & 
Small  (1963),  the  ozone  contribution  to  the  surface  pressure 
is  increased  by  about  58  per  cent  and  a  total  amplitude  is 
obtained  for  July  of  155  instead  of  the  present  113  jib. 

The  change  to  a  more  realistic  basic  temperature  profile  is 
therefore  the  main  factor  leading  to  the  present  reduced 
values . 

We  need  to  note  that  the  present  calculations  have  been 
based  on  classical  tidal  theory  which  ignores  mean  winds 
and  latitudinal  temperature  gradients.  Walterscheid  et  al. 
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(1980)  evaluated  semidiurnal  modes  under  both  classical  and 
non-classical  assumptions  and  found  that  the  introduction  of 
mean  winds  and  temperature  gradients  increased  the  solstitial 
values  of  the  (2,2,2)  surface  pressure  by  18  per  cent.  The 
vertical  structure  of  the  (3*3*4-)  mode  is  very  similar  to 
that  of  the  (2,2,2)  as  their  equivalent  depths  are  nearly 
equal,  being  7.66  and  7*85  km  respectively,  and  the  (3*3*4) 
surface  pressure  amplitude  is  likely  to  be  modified  by  a 
similar  order  of  magnitude.  Without  detailed  calculations, 
we  can  only  note  that  the  64  jxb  discrepancy  between  observation 
and  classical  theory  will  be  affected  by  non-classical  effects, 
possibly  by  as  much  as  20  pb. 

7.  Conclusions 

The  very  close  agreement  reported  (butler  &  Small,  1963) 
between  the  observed  and  calculated  amplitudes  of  the 
(3,3,4)  mode  is  found  to  be  fortuitous  as  it  depends  on  the 
choice  of  a  basic  temperature  profile  which  is  no  longer 
acceptable.  The  present  calculations  underpredict  the 
observed  solstitial  (3,3,4)  amplitude  by  about  64  pb,  i.e. 

36  per  cent.  The  significance  of  this  difference  is  uncertain 
as  observational  accuracy  has  not  been  defined,  but  previous 
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evaluation  of  diurnal  and  semidiurnal  components  (Groves, 
1981)  has  underpredicted  observed  values  and  it  is  therefore 
not  unreasonable  for  the  terdiurnal  component  to  be 
underpredicted  by  the  same  analytical  procedure. 
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